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SU (3) gauge theory with dynamical overlap fermions in the 2-index symmetric (sextet) represen- 
tation is considered. This model may be a viable model of the electroweak symmetry breaking 
sector along the lines of the walking technicolor paradigm. The number of fermion species is cho- 
sen such that the theory is expected to be below the conformal window. We will discuss how the 
£-regime and random matrix theory can be used to test whether at any given set of parameters {Nc, 
Nf, representation) the theory is in the conformal phase or indeed just below it. Quenched Monte 
Carlo results are included in the fundamental representation and also preliminary dynamical ones 
in the 2-index symmetric representation. 
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1. Introduction 

Model building of a strongly interacting electroweak sector, with or without a Higgs resonance, 
requires the knowledge of the phase diagram of non-abelian gauge theories for varying number of 
colors Nc, number of fermion flavors Nf, and representation R. For fixed Nc and R the theory is 
generically in the chiral symmetry broken phase for low Nf and the conformal phase for high Nf 
as long as asymptotic freedom is maintained, i.e. Nf is not too high. Certain models require Nf to 
be just below the conformal window along the lines of the walking technicolor paradigm [|l|] and 
the knowledge of the critical Nf separating the two phases is essential. 

Mapping out the phase diagram in the space of A'^^., Nf and R is an interesting problem on its 
own and can be useful for model builders with different motivations such as unparticles. We are 
first and foremost concerned with the Higgs mechanism though. 

In this context the parameters Nc, Nf and R are not only restricted by the phase diagram but also 
by electroweak precision data and the symmetry breaking pattern necessary for generating masses 
for the W and Z bosons. Consistency with electroweak precision data requires a small S'-parameter 
while the simplest symmetry breaking pattern is the one which generates exactly 3 Goldstone 
bosons with no (techni)pions left over after the massive gauge bosons acquired their masses. The 
simplest model fulfilling these requirements is SU (3) gauge group with Nf = 2 fermions in the 
2-index symmetric (25) representation which is the topic of our study. 

In a numerical simulation at finite volume, finite lattice spacing and (usually) finite quark 
masses it is a non-trivial task to determine whether the theory is conformal in the continuum, 
massless quark and infinite volume limits, or chiral symmetry is broken just as in QCD. In section 
^ a method is introduced that is capable of distinguishing the two phases based on the behavior of 
the low-lying Dirac eigenvalues. If the theory is QCD-like chiral symmetry is spontaneously broken 
and random matrix theory (RMT) will predict the distribution of low-lying eigenvalues in the e- 
regime whereas in the conformal phase chiral symmetry is unbroken and the spectral density 
of the Dirac operator goes to zero around A = 0. One particular advantage of RMT is that it works 
for finite (but small of course) quark mass. This method of distinguishing phases with the help of 
the Dirac spectrum has been applied for dynamical staggered fermions in [^] which complements 
Schrodinger functional and finite temperature based investigations of similar dynamical staggered 
models 

2. Perturbative expectations 

As is well-known the 2-loop perturbative jS -function ^ can be used to estimate the critical 
Nf value above which the theory is conformal and below which chiral symmetry is spontaneously 
broken [^. If the first two coefficients are j3i and ^2, asymptotic freedom requires j3i < otherwise 
the theory is free in the continuum. If j8i < and j32 > a non-trivial zero of the j3 -function exists 
hinting at a non-trivial IR fixed point. However if the fixed point value of the coupling is too large 
chiral symmetry is spontaneously broken before the flow in the IR can reach the would-be fixed 
point. It is nevertheless expected that a critical NJ" value exists above which the theory is really 
conformal even non-perturbatively. 
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The value NJ" can be estimated in the ladder approximation by the requirement that the 
anomalous dimension of y/i//^ reaches 7 = 1 ||8|]. Using this bound, the conformal window for 
SU (2) and representations j = 1 /2, 1 and 3/2 is expected to be 8 < Nf < 1 1, none (for integer Nf) 
and Nf = I respectively, and no window for j > 3/2. For SU (3) and fundamental, adjoint = 2A 
and IS representations the conformal window is expected to be 12 < Nf < 16, none (for integer 
Nf) and Nf = 3 respectively. 



3. Dirac spectrum 

To what extent the perturbative expectations of the previous section are justified is an open 
question in general. Non-perturbative tests of these expectations have been performed for various 



gauge groups, flavor number and representations f^, g, ^, |10|, [12|, |13|| using various methods. 

The low-lying spectrum of the Dirac operator is sensitive to the IR dynamics of the theory and 
shows characteristically different behavior in the conformal and QCD-like phases. Its measurement 
is in principle straightforward in a lattice simulation hence it is a good candidate to distinguish the 
two phases. 

3.1 Chirally broken phase, £-regime, random matrix theory 

If chiral symmetry is spontaneously broken, the Banks-Casher relation connects the spectral 



density p(A) of the Dirac operator around zero to the chiral condensate []14[], 



r=limlim^. (3.1) 

It also implies that the low-lying eigenvalues are dense in the sense that the average spacing is 
inversely proportional to the volume, 

AA = ^ . (3.2) 

It has been suggested long ago that if the bare parameters jS, m are tuned to the £-regime, i.e. 
such that rrij^ < < fn the low-lying Dirac spectrum follows the predictions of a random matrix 
theory [15, |l^]. The corresponding random matrix model is only sensitive to the pattern of chiral 



symmetry breaking, the topological charge and the rescaled fermion mass once the eigenvalues are 
also rescaled by the same factor ZV. 

More precisely, random matrix theory provides analytic formulae for the microscopic spectral 
density 

P5(C) = ^p(j7)=Iw(C), (3.3) 

and the individual eigenvalue distributions Pk{Q where ^ = AZV. The distributions Pk{Q only 
depend on pL = niLV, Nf and the topological charge v. The value of £ can be obtained by using jj. 
as the fitting parameter to have 

(0) = iM , (3.4) 
u m 
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where the left hand side is calculated in random matrix theory at a fixed charge v while the right 
hand side is measured in the simulation in the given sector v. Which eigenvalue Xu and which 
sector V is used is arbitrary in principle (as long as k is not too large, say, ^ = 1,2,3) and the quality 
of the whole procedure may be characterized by the (in)consistency of the obtained £ = iJ./{mV) 
values for various k and/or v. 

A more stringent test is the comparison of pk{Q between the random matrix theory predictions 
and the simulation once a consistent £ and corresponding jj. have been obtained from the above 
fitting procedure. The agreement is only expected for the first few eigenvalues. 

In order to see the effects of dynamical quarks the lowest eigenvalue should be larger than the 
fermion mass m. Otherwise the simulation is effectively quenched and random matrix theory will 
only agree at Nf = 0. 

Out of the two requirements of the £-regime, nijiL < 1 can be satisfied by tuning the fermion 
mass to a small value at any L. However the second requirement, fj^L > 1 , is largely independent 
of m provided it is small enough and puts a lower bound on L. As the lower edge of the conformal 
window is approached from below, /;;: is expected to decrease and eventually will vanish as the 
theory becomes conformal. Hence the fj^L > 1 condition will be more and more difficult to satisfy 
and larger and larger lattices will be needed the closer the theory is to the conformal window. As 
a result the study of nearly conformal (or walking) technicolor models is very challenging in the 
e-regime. 

It should be noted that the requirement f^^L > 1 is valid up to numerical constants only. From 
the behavior of the rotator and Goldstone spectrum of the chiral Lagrangian it can be made more 
precise as fj^L > 1/ \/27r which is the requirement of these two spectra to separate from each other. 
In fact we will see that in some cases RMT gives a good description even if fj^L < 1 which is 
probably due to the above numerical constant 1 / Vln = 0.3989... being smaller than 1. 

3.2 Conformal phase 

In the conformal phase no scale is generated and £ = 0. The spectral density of the Dirac 
operator around A ~ behaves as 



for massless quarks in the continuum and at infinite volume. Here 7 is the anomalous dimension of 

The exact dependence of 7 on the conformal fixed point coupling is in principle calcula- 
ble in perturbation theory since is presumably not large (otherwise chiral symmetry would be 
spontaneously broken). Certainly, 7 ~ gl. Of course only a non-perturbative treatment can decide 
whether there is room for a fixed point coupling which is large enough to be significantly different 
from perturbation theory and small enough so that chiral symmetry is not broken. 

How the ( p^ ) behavior is modified by finite volume and finite quark mass is an open question 
that we hope to address in the future. Certainly, in the free g^ = case the average eigenvalue 
spacing is inversely proportional to the linear size L of the box. The characteristic feature that is 
expected to hold even for a finite g^, > is that the average eigenvalue spacing for small eigenvalues 
is much less dense than in the chirally broken case where it is inversely proportional to the 4- volume 



(3.5) 



V; see (3.2). 



4 



Nearly conformal electroweak sector with chiral fermions 



Daniel Nogradi 



If the volume is too small, the chiral condensate is squeezed out of the box and the theory 
behaves perturbatively even in the case when chiral symmetry is broken in an infinite volume. 
Hence great care is needed not to confuse a small volume chiral symmetry breaking and a (small 
or large volume) conformal theory which is also behaving more-or-less perturbatively. 



4. Our model, SU{3) with A^^ = 2 in 25 representation 

The simplest example of a model that - according to the perturbative expectations - is just 
below the conformal window, has a relatively low Nf value so that the S'-parameter is relatively 
small, and has precisely 3 Goldstone bosons is Nc = 3, Nf = 2 and R = 2S. This model has been 
studied in [[u], |l2[ ll] using Wilson fermions on rather small lattices and it was found to be already 
in the conformal window although it was indicated that more complicated possibilites are also 
allowed by the data. 

Since exact chiral symmetry is important both for QCD-like and conformal theories we chose 



to use overlap fermions [17]. The simulation has to be carried out at a fixed topological charge. 



There are two methods available for simulating dynamical overlap fermions at fixed topology. 



One is the reflection/refraction algorithm [18] but always reflecting on the topological boundary. 



The other is employing a pair of extra Wilson fermions to suppress exact zero modes thereby 



suppressing tunnelling between sectors [19]. We used the second method, which is much faster, in 
this study. 



5. Preliminary results 

5.1 Quenched simulations 

In order to see how well actual simulations agree or disagree with the predictions of random 
matrix theory (RMT) we have tested the RMT predictions for Nf = and overlap valence quarks in 
the fundamental representation since in the quenched approximation chiral symmetry is guaranteed 
to be broken. This setup is identical to [ pO| ] but actual eigenvalue distributions were not presented 
there. Since valuable information can be gained from these we decided to redo this analysis on 12^ 
lattices at j8 = 5.8458 which corresponds to a lattice size of L = 1.49 fm. All our parameters were 
the same as in [20]. In particular we used the Wilson gauge action and an unsmeared overlap oper- 
ator. Our ensemble consists of 1500 configurations. The results of [^] for expectation value ratios 
(A,) / (Ay) have been reproduced within 1-sigma precision with occasional 1.2-sigma deviations. 

The distribution of the ^th eigenvalue in various topological sectors \Q\ are shown in figure |l] 
together with the RMT predictions once £ has been fitted from k = I and 2 = as described by 



(3.4). Clearly, not only the expectation values (A,) follow the RMT predictions but also their width. 



This level of agreement with RMT was not expected since fj^L < 1 for this ensemble, however as 



indicated in section |3jJ a more accurate requirement of the £-regime is fjiL > Xj^/lM which this 
ensemble does fulfill. 

A similar comparison for the 2-index symmetric representation in the quenched approximation 
is ongoing. Since RMT is only sensitive to the topological charge, Nf, jJ. and the pattern of chiral 
symmetry breaking, the same random matrix model is expected to describe this representation as 
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Q = 0, k.l = 1, k.1 Q = 2, k.1 
0.6 I , 1 1 , 1 1 1 0.45 I 1 1 , 1 , 1 , 1 , 1 0.4 I , 1 , 




01234567 0123456789 10 2 4 6 8 10 12 

Q = 0, k = 2 Q = 1,k.2 Q = 2, k.2 




0123456789 10 11 02468 10 12 14 02468 10 12 14 



Figure 1: Rescaled quenched eigenvalue distributions pk{Q in the fundamental representation and Nf ^0 
RMT predictions for |2| = 0, 1 , 2 and ;t = 1 , 2 



the one used for the fundamental. Agreement with the same RMT using a different representation 
than fundamental will be a non-trivial check of its universality. 

5.2 Dynamical simulations 

In the 2-index symmetric representation three dynamical ensembles were generated on 6^ lat- 
tices using the tree-level improved Symanzik gauge action at j8 = 4.850, 4.975 and 5.100 and 
Nf = 2 flavors of massive quarks with m = 0.05. The negative Wilson mass in the overlap operator 
was mw = — 1.3 and 2 levels of stout smearing with smearing parameter p =0.15 have been ap- 



phed. The topology change suppressing action of [19] was used with mass M = 0.2 for the ghost 
Wilson fermions and only the topological sector Q = was sampled. Since chiral symmetry is 
preserved by overlap fermions at finite lattice spacing Nf = 2 RMT is applicable. Fitting £ from 
the average first eigenvalue as in ( |3.4| ) one obtains 0.083(4), 0.084(4) and 0.080(4) in lattice units 
for the three j8 values respectively. 

The eigenvalue distributions look qualitatively the same for the three j8 values and the first 3 
eigenvalues are plotted in figure || for j8 = 4.850 together with the Nf = 2 RMT predictions after 
rescaling both the eigenvalues and the quark mass. 

Clearly, the RMT predictions are very far from the simulation results, neither the averages nor 
the widths follow the RMT curves. This may be due to several reasons the most likely of which is 
small volume. We have not measured either niji or fji so it is not clear if the simulation was in the 
£-regime at all. Simulations on larger volumes as well as measurements of niji and fji are ongoing. 

If the larger volume simulations agree with the above conclusion the 2-index symmetric rep- 
resentation for gauge group SU (3) and Nf = 2is already in the conformal window. 

6. Conclusions and outlook 

Needless to say that the results on the 2-index symmetric representation are preliminary. The 
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k=1 k.2 k.3 




Figure 2: Rescaled dynamical eigenvalue distributions in the 2-index symmetric representation and 

Nf = 2 RMT predictions for g = and /t = 1 , 2, 3 



quenched fundamental representation simulations shows that the RMT predictions are very precise 
for the first few eigenvalues once the volume is large enough. A similar conclusion is expected 
for the quenched 2S representation which will in addition test the universality of RMT. Presently, 
the result of the 6^ dynamical 2S representation simulation are preliminary and the deviation from 
RMT is thought to be due to small volume. Larger volume simulations are ongoing for both the 
quenched and fully dynamical cases. 
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